Impact of the matter density uncertainty on the dark energy reconstruction 
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In this paper we study the impact of the fractional matter density uncertainty in the reconstruction 
of the equation of state of dark energy. We consider both standard reconstruction methods, based 
on the dynamical effect that dark energy has on the expansion of the Universe, as well as non- 
standard methods, in which the evolution of the dark energy equation of state with redshift is 
inferred through the variation of fundamental couplings such as the fine structure constant, a, or 
the proton-to-electron mass ratio, ^. We show that the negative impact of the matter density 
uncertainty in the dark energy reconstruction using varying couplings may be very small compared 
to standard reconstruction methods. We also briefly discuss other fundamental questions which 
need to be answered before varying couplings can be successfully used to probe the nature of the 
dark energy. 
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I. INTRODUCTION 

Multiple observations in the last decade have provided 
very strong evidence for a recent acceleration of the uni- 
verse [U [11 El SISl IS] • In the context of Einstein General 
Relativity such an acceleration can only be explained by 
an exotic dark energy component violating the strong en- 
ergy condition. From a purely phenomenological point of 
view the cosmological constant appears to be the simplest 
dark energy candidate. Still, a dynamical scalar field is 
expected to be a more plausible explanation particularly 
in face of the very large discrepancy between the observa- 
tionally inferred vacuum energy density and theoretical 
expectations. 

If the dark energy is dynamical a fundamental question 
immediately arises regarding the characterization of the 
evolution of its properties with redshift and in particular 
of its equation of state. Standard methods to reconstruct 
the dark energy equation of state as a function of the red- 
shift rely on the dynamical effect that dark energy has 
on the expansion of the universe [7j. However, dark en- 
ergy is dynamically relevant mainly at recent times which 
makes the task of accurately determining the evolution 
of its equation of state at z > 1 an almost impossible 
one, at least using standard methods. 

On the other hand, cosmological variations of funda- 
mental couplings can be probed over a wide redshift 
range. At high redshift cosmic microwave background 
temperature and polarization anisotropics O El UHl [HI 
|12 , 13, 14, 15J and hght element abundances [9} [Ml [TBI [T7] 
constrain the value of a at z '--^ 10^'^ and z 10^ to be 
within a few percent of its present day value. Positive 
results suggesting a cosmological variation of the fine- 
structure constant, a, and the proton-to-electron mass 
ratio, /i, at about the 10~^ have been reported in the 
redshift range z = l- 4[THl[iniini[lIl- Unfortunately, 



other analysis have found no evidence for such varia- 
tions im [231 mi- At even lower redshifts laboratory 
experiments and the Oklo natural nuclear reactor pro- 
vide very stringent limits on the time- variation of a and 



If the dark energy is described by a quintessence field, 
(j), non-minimally coupled to the electromagnetic field 
[S[30l[ai[321[33[3lllMl[Mll37l[MllMl[iniS^ then 
the dynamics of a is coupled to the evolution of (j). It 
was shown 42J that, under certain assumptions, varying 
couplings may be used to probe the nature of dark energy 
over a larger redshift range than that spanned by stan- 
dard methods (such as supernovae [H [H |2l Hj or weak 
lensing (131 HH US]). However, a perfect knowledge of 
ri.,„o was assumed which can not be realized in practice. 
In the present paper we eliminate this shortcoming by 
studying the impact that the matter density uncertainty 
has on the dark energy reconstruction using varying cou- 
plings. 



This paper is organized as follows. In Section II we 
start by quantifying the impact of the fractional matter 
density uncertainty in the determination of the evolution 
of the dark energy equation of state with redshift using 
standard reconstruction methods. In Section III we per- 
form a quantitative analysis of the same problem now 
assuming that varying couplings are used to probe the 
nature of dark energy. We consider a broad class of mod- 
els for the evolution of a and fi where the gauge kinetic 
function is a linear function of a quintessence-type real 
scalar field described by a Lagrangian with a standard 
kinetic term and a scalar field potential, V {(/)). We con- 
clude in Section IV with a brief summary of our results 
and a discussion of future prospects. 
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Throughout this paper we shall use units with h = c = 
SttG = 1 and a metric signature (-1-, — , — , — ). 
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II. STANDARD DARK ENERGY 
RECONSTRUCTION 

Consider a flat homogeneous and isotropic Friedmann- 
Robertson- Walker universe whose dynamics is described 

by 



HH' 



tt2_P 
\{P + P) , 



(1) 

(2) 



where p is the total density, p is the total pressure, a 
is the scale factor, a prime represents a derivative with 
respect to In a, H — a /a is the Hubble parameter and a 
dot represents a derivative with respect to physical time, 
t. Eqs. ([ij and ^ can also be combined to give 

Eq. ([3| tell us that a precise knowledge of the universe 
dynamics (and therefore of the evolution of H) is all that 
is required in order to determine the evolution of w. If the 
universe contains only minimally coupled fluids energy- 
momentum conservation implies that the equation 



p' + 3(p + p) = 0, 



(4) 



is satisfied not only by the mixture but also by each in- 
dividual fluid. In this paper we will make the simpli- 
fying assumption that the universe is constituted solely 
by matter and dark energy with energy densities pm and 
/50, thus neglecting the residual (at recent times) radia- 
tion component. Consequently, p = pm + P4, and p = Pcf, 
(taking pm = 0) so that 



1 -t- n,n/n^ 



W^{1 - flrn) , 



(5) 



where flm = Pm/p, ^ci> — P4>/ P and Q,m + ^4, — ^ since 
we are considering a flat universe. 

Standard methods to reconstruct the dark energy 
equation of state rely on the dynamical effect that dark 
energy has on the expansion of the universe. They are 
faced with two main limitations. The first is related to 
the fact that the dynamics of the universe is not known 
with infinite precision. Hence, at each redshift, z, the 
value of w will have an uncertainty Aw that will trans- 
late into a much larger uncertainty in Aw^ 



Aw0 



1 



Aw 



(6) 



For the moment we are assuming no uncertainty in f2m 
(i.e. = 0). 

The ratio grows rapidly with redshift at recent 

times and consequently we expect that Aw^ ^ Aw for 
z > 1 (note that 1 -I- z = 1/a). If we assume that the 
dark energy is a cosmological constant then 

il„j ^ il,„ ^ _^ ^ i7„,o(l + 



1 



AO 



1 + 0.37(1 -I- z) = 



(7) 
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FIG. 1: Evolution with redshift, z, of the uncertainty in the 
value of assuming a background ACDM model with Q.mO = 
0.27, Q,^o = SIao = 0.73 and an observational uncertainty 
Ailmo = 0.01 in the fractional matter density. The dashed 
line represents the correct value of the dark energy equation 
of state parameter (equal to —1 independently of z, in this 
particular case). 



which becomes a very large factor at high redshift (here 
we took ri„jo = 0.27 and ilAo = 0.73 as favored by the 
five-year WMAP results [6,). The subscript '0' indicates 
that the variables are to be evaluated at the present 
epoch. 

A second source of uncertainty in is due to the fact 
that we do not know ri„jo with infinite precision [1^1 SZl 
I48j . This uncertainty would be present even if we had 
a perfect knowledge of the evolution of w (in which case 
Aw would be zero). Although this ideal situation cannot 
be achieved in practice, we will consider it for the sake 
of illustration. Still, a perfect knowledge of the dynamics 
of w would not imply that the evolution of w^ could also 
be determined with arbitrary precision. An uncertainty 
in the value of matter density, Ap^, will be reflected in a 
corresponding uncertainty in the value of p^, since we are 
unable to determine, with absolute certainty, whether or 
not a fraction of the matter density is really dark matter 
or dark energy. The corresponding uncertainty in the 
dark energy equation of state parameter is equal to 



Aw0 



_HAf^m_ 



Ap„ 



(8) 



where Afim = ^Pm/ p > and the approximation is valid 
only if Aw0 If the dark energy is a cosmological 

constant then we have 



nr. 



(9) 



with f2 



f^mO + f^Ao(l + 2^) 

1 — f2„, . On the other hand we assume 
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that 



(10) 



In Fig. 1 we plot the evolution with redshift, z, of 
the uncertainty in the dark energy equation of state pa- 
rameter, [w^ — Aw^,w^ + Awc/,], assuming again that 
^mo = 0.27, flcfio = f^Ao = 0.73 and an observational 
uncertainty Aflmo = 0.01 in the fractional matter den- 
sity. The dashed line represents the correct value of the 
dark energy equation of state parameter (equal to —1 
independently of z, in this particular case). The figure 
clearly shows a small uncertainty in the value of 
z = originates a very large uncertainty in the evolu- 
tion of the dark energy equation state, specially at high 
redshift. 



III. DARK ENERGY RECONSTRUCTION 
USING VARYING COUPLINGS 

In this section we shall consider a class of models mod- 
els where a quintessence field is non-minimally coupled 
to the electromagnetic field. These models are defined 
by the action 



S = 



where 



C — + Ccf,F + -Mother , 



V{(t>) is the scalar field potential. 



1 



C4,F = --BF{4>)F^,,F^''' . 



(11) 



(12) 



(13) 



(14) 



Bp^cj)) is the gauge kinetic function, F^^ are the compo- 
nents of the electromagnetic field tensor and £othcr is the 
Lagrangian density of the other fields. The fine-structure 
constant is given by 



ao 



Bf{<P) 



(15) 



and, at the present day, one has Bf{4>q) = 1. 

We will also assume that the gauge kinetic function is 
a linear function of 6 so that one has 



da 
a 



(16) 



where 5a = — a, Scj) ~ — (j) and C is a constant. 

The evolution of ((> induced solely by its coupling to 
electromagnetically interacting matter is so small (given 



Weak Equivalence Principle constraints [301 HSl ISO] ) that 
the resulting time variation of a can be neglected. Hence, 
throughout this paper we shall assume that the dynamics 
of ip is fully driven by the scalar field potential, V{(j)) 
(and damped by the expansion). Throughout this paper 
we shall also neglect the spatial variations of cj) which 
was shown to be a good approximation in this context 

gniisii 52 . 

The relation between the variations of a and /i is model 
dependent but, in general, we expect that 



/i a 



(17) 



where i? is a model dependent constant (see [53J [SH [S5J 
|56l [57l [58l |59] for detailed discussion of specific models). 

Taking into account that p = pm + P4,, P4, = 0^/2 
F(0), = </.V2 - y(0), p„, = and ^ = <j)'H, Eqs. |l 
and ^ can be written as 



(Pm + P<t>) 



-3iJ 



or alternatively 







H 



"0 



where a = P(ji/{pa^mo) so that 



Substituting Eq. ( 20 1 into Eq. (|21j) one obtains 

2 



co + 1 



"0 



(18) 
(19) 

(20) 
(21) 

(22) 

I 

(23) 



Note that the linearity assumption given by Eq. (16 1 
implies that (f)' l4>'o = a' /a^. The dark energy equation 
of state parameter is given by 



= -1 



P<t>+P<t> 



= -1 - 



and consequently 



1) 



-3a-o 



= -1 - 



1 - SI 



Eq. ( 24 1 can also be written as 

W0 = -IH = -1 



mO 



3f}w 



which implies that 



1 



00 



a 



a; 



(24) 



(25) 



(26) 



(27) 
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FIG. 2: The evolution with redshift, z, of the equation of 
state parameter, u/^, determined using Eqs. (27 1 and (281 
and taking wo ~ —0.68, Jlmo = 0.27 + &VLmO with iVtmn = 
-0.04,-0.02,-0.01,0.01,0.02,0.04 (n<^o = 1 - f^mo). The 
blue stripe represents the interval — Aw^ , + Ait;^] with 
the lower and upper limits of the interval calculated using 
f^mo = 0.26 and SlmO = 0.28 respectively. Here we consider 
a model with constant a' (upper panel) and another with 



.1/2 



(lower panel). 



where 



4> 



P 



(28) 



In this paper our aim is to isolate the effect of the un- 
certainty in ri„io on our ability to reconstruct the evolu- 
tion with redshift of the dark energy equation of state 
and consequently we shall take the ideal situation in 
which the evolution of a with redshift is known with in- 
finite precision (note that forecasts taking into account 
expected measurement uncertainties for the next gener- 
ation of spectrographs have already been presented else- 
where IH]). Hence, here we assume that the uncertainty 
in /^(z) is completely due to the uncertainty in ^^{z) 
(or equivalently in <7{z)). 

In Fig. 2 we plot the evolution with redshift, z, of the 
equation of state parameter, Wrf,, determined using Eqs. 
( [Zf| and ([28| and taking wq = -0.68, = 0.27 -I- 

d^mo with S^mQ = -0.04,-0.02,-0.01,0.01,0.02,0.04 



(fi^o = 1 — f^mo)- The blue stripe represents the interval 
[w^ — Aw0, + Aw^] with the lower and upper limits of 
the interval calculated using O^o = 0.26 and ri,,„o = 0.28 
respectively. For the sake of illustration we consider a 
model with constant a' (upper panel) and another with 
a' cx a}/'^ (lower panel). However, we have verified that 
our main conclusions will not depend on these particular 
choices for the evolution of a. 

We clearly see in Fig. 2 that the error bars on the 
dark energy equation of state reconstruction (blue stripe) 
due to the uncertainty in the value of f^mo do not tend 
increase with redshift. We can see that is indeed the 
case by considering a particular set of models in which a 
becomes considerably greater than ctq at large z. When 
z > 1 we have aa? < 1 and consequently 



aa 



where 



6 



Co a 



-3(1 -I- W0o) 



The uncertainty in the value of ^0 is given by 

(1 + W0o)^n^ = Aw0o • 
l4o| 



(29) 
(30) 

(31) 
(32) 



If the value of a at high redshift is mainly determined by 
its evolution at z > 1 then Eqs. (27 1, (29) and (30 1 imply 



that is roughly proportional to and consequently 
dftp/fj^jy — d^o/Co for z > 1. On the other hand, from 
Eq. ([27| we have that dw^ = fcpidw4o + {l + w^o)df4,/f^) 
and we may now estimate the high redshift uncertainty 
in the equation of state parameter in these models to be 



Aw0 Aw 



ACo 



0. 



(33) 



Although this approximation is not valid at low red- 
shift or while cr ^ uo it clearly shows that, contrary to 
standard reconstruction methods, the uncertainty in the 
equation of state parameter due to the uncertainty in 
flmo does not blow up at high redshift (assuming that 
the lower limit of the interval [w^o — Aw^o, w^o + Aw^o] 
is significantly above —1). This is hardly surprising since 
if we were able to track the evolution of on a with 
redshift with arbitrary precision (up to a normalization 
factor) then we would also have a perfect of knowledge 
of w^{z). 

However, if we consider values of w^q close to —1 then 
the distance between adjacent lines in Fig. 2 (upper and 
lower panels) tends to increase with redshift. This is 
related to the fact that as we approach w^o = — 1 we 
get closer to the cosmological constant case for which 
— —1 independently of z. Still, we note that an 
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observation of a variation of a at z = immediately 
leads to an upper bound on w^o [50] . In particular, from 
Eqs. ((T6|) and p6|) one obtains 



w^o = -1 



(34) 



where C ^ 1.8 x 10~* considering a variation of a alone 
(i.e. fixing all other gauge couplings). The constrains 
on C are in general much stronger when more realistic 
models, based on the unifications of gauge couplings, are 
considered [50]. Hence, if ^ then w^o cannot be ar- 
bitrarily close to —1. The relation between experimental 
limits at z = and varying coupling constraints higher 
redshifts is of course model dependent. However, if 4> 
is always rolling down a monotonic potential then low- 
redshift constrains on the evolution of a may already 
provide stringent limits on the corresponding variations 
at high-redshift [IT. 

The upper panel of Fig. 2 shows that the error bars 
due to the uncertainty in the value of flmo can become 
significantly smaller at large redshift. Note that Eq. Q 
implies that for a constant 0' the energy density of the 
scalar field is expected to track the background density 
at early times since 



oc p . 



(35) 



If l^gl is not too small then ^ deep in the matter 
era irrespectively of the initial conditions. This is the 
main reason why the uncertainty becomes very small at 
large redshift in the upper panel of Fig. 2. We have also 
verified that if </>' (x a^^"^ (lower panel) then 
asymptotically at large z. However, in this case the con- 
vergence is slower than in the previous one. One should 
however bear in mind that future spectroscopic determi- 
nations are unlikely to be able to probe the dynamics of 
dark energy beyond z = 5 — 6. 



methods. Standard methods are very sensitive to this 
uncertainty since the ratio rim/fi^ is expected to be- 
come very large at high redshift. On the other hand, 
we have shown that the negative impact of the matter 
density uncertainty may be much smaller in the case of 
non-standard reconstruction methods in which the evolu- 
tion of the dark energy equation of state with redshift is 
inferred through the variation of fundamental couplings 
(such as a or /i). 

There are however a number of other important ques- 
tions which still need to be answered positively if varying 
couplings are to be successfully used to infer the dynam- 
ics of dark energy: I. Do a (/i) evolve with cosmic time 
and are such variations observable ? H. Is dark energy 
a quintessence scalar field, (j), described by a Lagrangian 
with a standard kinetic term and an effective potential 
which is a function of (j) alone ? III. Are the variations of 
a (/i) linearly coupled to ? 

Presently, the answer to the first question is still con- 
troversial. Although a number of results suggest a cosmo- 
logical variation of a and p, in the redshift range z = 1 — 4 
other analysis have found no evidence for such variations. 
However, this situation is expected to be fully clarified 
in the next few years with the next generation of high 
resolution spectrographs ^U\. Unfortunately, questions 2 
and 3 may be easier to answer negatively than positively 
(for example if standard and non-standard reconstruction 
methods give incompatible results). Also, note that the 
varying-couplings method is in general affected by the 
observational uncertainty in the present-day equation of 
state wq, which is necessary to calibrate the relation be- 
tween (j)' and a' |36j . Still, given the strong limitations 
of standard reconstruction methods alternative promis- 
ing approaches, such as those based on the time- variation 
of fundamental couplings, deserve to be further investi- 
gated. 



IV. CONCLUSIONS 

In this paper we quantified the impact of the frac- 
tional matter density uncertainty, AflmOi in the recon- 
struction of the equation of state of dark energy con- 
sidering both standard and non-standard reconstruction 
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